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Coenzyme B12 (AdoCbl) is one of the most biologically active forms of vitamin B12, 
and continues to be a topic of active research interest. The mechanism of Co-C bond 
cleavage in AdoCbl, and the corresponding enzymatic reactions are however, not well 
understood at the molecular level. In this work, time-dependent density functional theory 
(TD-DFT) has been applied to investigate the photodissociation of coenzyme B12. To 
reduce computational cost, while retaining the major spectroscopic features of AdoCbl, 
a truncated model based on ribosylcobalamin (RibCbl) was used to simulate Co-C 
photodissociation. Equilibrium geometries of RibCbl were obtained by optimization at the 
DFT/BP86/TZVP level of theory, and low-lying excited states were calculated by TD-DFT 
using the same functional and basis set. The calculated singlet states, and absorption 
spectra were simulated in both the gas phase, and water, using the polarizable continuum 
model (PCM). Both spectra were in reasonable agreement with experimental data, and 
potential energy curves based on vertical excitations were plotted to explore the nature 
of Co-C bond dissociation. It was found that a repulsive 3 (aco-c -> ^co-c' triplet state 
became dissociative at large Co-C bond distance, similar to a previous observation for 
methylcobalamin (MeCbl). Furthermore, potential energy surfaces (PESs) obtained as a 
function of both Co-CRjb and Co-N| m distances, identify the S^ state as a key intermediate 
generated during photoexcitation of RibCbl, attributed to a mixture of a metal-to-ligand 
charge transfer (MLCT) and a a bonding-ligand charge transfer (SBLCT) states. 

Keywords: coenzyme Bi 2 , Co-C bond, photodissociation, ribosylcobalamin, time-dependent density functional 
theory 



INTRODUCTION 

Vitamin B12 derivatives (Figure 1) are a group of organometal- 
lic complexes that act as biologically active cofactors in many 
enzymatic reactions (Dolphin et al, 1982; Banerjee, 1997, 1999, 
2001, 2003; Ludwig and Matthews, 1997; Krautler et al, 1998; 
Marzilli, 1999; Toraya, 2000; Matthews, 2001; Banerjee and 
Ragsdale, 2003; Toraya, 2003; Brown, 2005; Randaccio et al., 
2006, 2007). In addition to their pivotal roles in a variety of 
enzymatic processes, the B12 derivatives possess complex photo- 
physical and photochemical properties (Endicott and Ferraudi, 
1977; Endicott and Netzel, 1979; Rao and Symons, 1982; Chen 
and Chance, 1990, 1993; Sakaguchi et al, 1990; Chagovetz 
and Grissom, 1993; Grissom and Chagovetz, 1993; Lott et al., 
1995; Natarajan and Grissom, 1996; Kruppa et al, 1997; Walker 
et al, 1998a,b; Shiang et al, 1999, 2006; Yoder et al, 2001; 
Cole et al, 2002; Sension et al., 2004, 2005a,b; Harris et al., 
2007). The relatively weak organometallic Co-R bond (R = 
alkyl, hydroxyl, water, cyanide) in these compounds can undergo 
photodissociation under conditions of simple photon excitation, 
depending on the nature of the upper axial ligand. Cobalamin 
complexes with the alkyl axial ligands, such as enzymatically 
competent methylcobalamin (MeCbl) and adenosylcobalamin 
(AdoCbl), and their analogs, ethylcobalamin (EtCbl) and propy- 
lcobalamin (PropCbl), undergo photodissociation upon expo- 
sure to light (Walker et al, 1998a,b; Shiang et al, 1999; Yoder 



et al., 2001; Cole et al, 2002; Sension et al, 2004, 2005a,b; 
Harris et al, 2007). In certain cases, the nature of events fol- 
lowing photon excitation is additionally influenced by a wave- 
length of exciting light, as found for MeCbl (Harris et al., 
2007). On the other hand, the Co-R bond in cyannocobal- 
imin (CNCbl) and non-alkylcobalamin analogs, such as azido- 
cobalamin (N3CM) and aquocobalamin (H20Cbl + ), are rather 
photostable (Shiang et al., 2006). 

The photochemistry of cobalt corrinoids has been probed 
using various experimental techniques such as laser flash pho- 
tolysis (Endicott and Ferraudi, 1977; Endicott and Netzel, 1979; 
Chen and Chance, 1990; Chagovetz and Grissom, 1993; Lott 
et al., 1995), continuous wave (CW) photolysis (Chen and 
Chance, 1993), kinetic magnetic field effect (MFE) (Grissom 
and Chagovetz, 1993; Natarajan and Grissom, 1996), chemi- 
cally induced dynamic electron polarization (CIDEP) (Rao and 
Symons, 1982; Sakaguchi et al, 1990) as well as nuclear polar- 
ization (CIDNP) (Kruppa et al, 1997). Systematic studies of 
photolysis in B12 alkyl- derivatives (MeCbl, EtCbl, PropCbl, and 
AdoCbl) using pump-probe transient absorption spectroscopy 
found that photoproduct yields for alkylcobalamins are differ- 
ent when excited at 400, or 520/530 nm (Walker et al, 1998a,b; 
Shiang et al, 1999; Yoder et al, 2001; Cole et al, 2002; Sension 
et al, 2004, 2005a,b; Harris et al, 2007). MeCbl for example, pro- 
duces a mixture of 27% cob(II)alamin and 73% cob(III)alamin 
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FIGURE 1 | Upper: Molecular structure of vitamin B-12 derivatives where R 
= Me, Ado, Et, Prop, CN, OH or N 3 where Ri =CH 2 CONH 2 , R2 = 
CH 2 CH 2 CONH 2 , and R 3 = (C^bCONHC^CHIChyOPO^ . Lower: 
Structural model of RibCbl employed in present work (Rib refers to ribosyl 
with 5-hydroxyl group substituted by H atom). Reprinted (adapted) with 
permission from Jaworska et al. (2007). Copyright (2014) American 
Chemical Society. 



when excited at 400 nm as a result of partitioning between Co-C 
homolysis and heterolytic cleavage. Heterolytic bond cleavage 
products were identified as a metal-to-ligand charge transfer 
(MLCT) state upon excitation at 400 nm (Cole et al, 2002), 
and only metastable cob(III)alamin for 520 nm (Shiang et al., 
1999). In contrast, photoexcitation of AdoCbl produces wave- 
length independent photoproducts at 400 and 520 nm, that 
correspond to an intermediate state observed as decay over a 
14 ps time scale, of two relaxed radical pairs. These radical pairs 
were characterized as cob(II)alamin and an adenosyl-based rad- 
ical. This observation corresponds to Co-C homolysis, or a 
trapped excited state with a weak Co-C bond (Shiang et al., 
1999). Furthermore, although the photoproducts of AdoCbl 
were found to be insensitive to excitation wavelength, environ- 
mental effects have been shown to influence its photophysics 
(Yoder et al, 2001). 

The photodissociation of Co-C in cobalamins has also been 
investigated computationally by means of time-dependent den- 
sity functional theory (TD-DFT) (Jaworska et al., 2007; Lodowski 
et al, 2009, 2011; Kumar and Kozlowski, 2012). In the case of 



MeCbl, the key metastable photoproduct present during the pho- 
tolysis process was identified as an Si state with MLCT character 
(Jaworska et al., 2007), in accordance with experiment. Further 
analysis of the electron density map of Si and So states indi- 
cated an additional contribution from 0 bond, revealing that 
a bonding-ligand charge transfer (SBLCT) character was mixed 
with MLCT (Lodowski et al, 2009). To explain the nature of 
the Co-C photo-scission in MeCbl, the presence of a repulsive 
3 (ctco-c — > °co-c) triplet state was proposed based on TD-DFT 
analysis of electronically excited states along the stretched Co-C 
bond. The same approach has also been applied to EtCbl, and 
the presence of a repulsive triplet state was used to explain the 
energetic differences between MeCbl and EtCbl, as well as why 
the photolysis mechanism is wavelength dependent in the case of 
MeCbl (Lodowski et al., 2009). 

A similar study was performed for vitamin B12 (CNCbl) to 
explain its photostability and the nature of its low-lying excited 
states (Lodowski et al., 2011). Potential energy curves for low- 
lying excited states along the Co-CN coordinate reveal that CNCbl 
has a repulsive triplet state, but is not dissociative. It was found 
that the CNCbl potential energy surface (PES) of the Si state, 
when represented as a function of axial distances, had two energy 
minima. The first minima was located above the So minimum 
as an excitation with mixed jt — > Tt*/MLCT/SBLCT character, 
similar to the alkylcobalamin Si state. The second minimum 
was found at longer Co-Ni m and Co-CN bond lengths, and 
was characterized as ligand-to-metal charge-transfer (LMCT) 
excitation. 

Interestingly, the Co-OH bond in hydroxycobalamin (HOCbl) 
is dissociative upon light exposure above 300 nm (Shell and 
Lawrence, 2011), although it is among the non-alkylcobalamins. 
Photolysis of HOCbl has also been studied by TD-DFT, by focus- 
ing on the Co-OH cleavage (Kumar and Kozlowski, 2012). TD- 
DFT results suggest that the photoactivity of HOCbl is mediated 
by the repulsive Hn+dco — >- crJ 0 _ OH ) singlet state whose energy 
drops with Co-OH bond stretching, to yield the cob(II)alamin 
and hydroxyl radical. 

Although spectroscopic techniques have been extensively 
applied to alkyl- and non-alkylcobalamins, the mechanism of 
Co-C bond dissociation is still not well understood. The aim of 
this study is to provide the further insight into the photodis- 
sociation mechanism of the Co-C bond in the coenzyme B12 
(AdoCbl, Figure 1) by TD-DFT computations. To accomplish 
this we employed the simplified structural model of AdoCbl, 
(RibCbl, Figure 1), in which the Ado group was simplified 
to a ribosyl (Rib) moiety. The intermediates involved in the 
photodissociation process were identified by systematically ana- 
lyzing the nature of low-lying excited states as well as their 
changes along stretched Co-Ciub coordinate. Finally, in order 
to provide a more accurate picture of key intermediates, the 
PESs of ground state and low-lying singlet states were ana- 
lyzed as a function of both Co-CRib and Co-Ni m axial bond 
distances. 

COMPUTATIONAL DETAILS 

All calculations reported in this work were based on non- 
local DFT with the gradient-corrected Becke-Perdew (BP86) 
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(Perdew, 1986; Becke, 1993) functional and the TZVP basis 
set, as implemented in the Gaussian 09 (Frisch et al., 2009) 
or TURBOMOLE (Ahlrichs et al., 1989; Treutler and Ahlrichs, 
1995; Furche and Ahlrichs, 2002; Furche and Rappoport, 
2005; TURBOMOLE : ) suites of programs for electronic struc- 
ture calculations. The BP86 functional has been previously 
selected as the most appropriate for predicting both the struc- 
tural and electronic properties of B12 cofactors, including Co- 
C bond dissociation energy (BDE) (Jensen and Ryde, 2003; 
Rovira et al., 2004; Kuta et al., 2006; Kozlowski et al., 2007; 
Rovira and Kozlowski, 2007; Galezowski et al, 2008; Lodowski 
et al, 2009, 2011; Kornobis et al, 2011, 2013). To account 
for environmental effects on geometries and electronic proper- 
ties, ground states of RibCbl model (Im-[Co m (corrin)]-Rib + , 
see Figure SI) were computed in the gas phase as well as in 
water via the Polarizable Continuum Model (PCM) (Miertus 
et al., 1981; Cammi and Tomasi, 1995) as implemented in 
Gaussian 09. 

Vertical excitations were obtained in both the gas phase and 
water solvent. Low-lying excited states were calculated from their 
corrdepsonding ground state geometries by TD-DFT, and both 
singlet and triplet excitations have been calculated to forty states 
(Tables S1-S4). To simulate the photodissociation process of 
RibCbl, the Co-Cmb bond was stretched with a step size of 
0.05 A and the DFT/BP86/TZVP optimization was repeated at 
each point, followed by TD-DFT calculations at the BP86/TZVP 
level of theory. In addition, PESs of the ground state and key 
singlet states were plotted by stretching axial bonds in step of 
0.05 A. At each grid point, the geometry was optimized by 
DFT/BP86/TZVP, and the vertical excitation was obtained for Si 
and S2 states. 

RESULTS AND DISCUSSION 
STRUCTURAL MODEL OF RibCbl 

The initial molecular structure of RibCbl was obtained from 
the high-resolution crystal structure of AdoCbl (Ouyang et al., 
2004). The full structure of AdoCbl (Figure 1) consists of Co 111 
equatorially cooridinated via four nitrogen atoms to the cor- 
rin ring. The lower axial ligand, dimethylbenzimidazole (DBI) 
coordinates to the central Co atom via N ax , with the other N 
atom bound to the side of the corrin plane by a nucleotide 
loop. In the actual AdoCbl structure, the upper Ado ligand is 
5'-deoxy-5'-adenosine. 

To reduce computational cost, calculations were performed 
using the simplified structural model of RibCbl, denoted 



1 TURBOMOLE has been designed by the QuantumChemistryGroup, 
University of Karlsruhe, Germany, since 1988. The following members of 
the group have made contributions: Reinhart Ahlrichs, Michael Bar, Hans- 
Peter Baron, Riidiger Bauernschmitt, Stephan Bocker, Nathan Crawford, Peter 
Deglmann, Michael Ehrig, Karin Eichkorn, Simon Elliott, Filipp Furche, 
Frank Haase, Marco Haser, Christof Hattig, Arnim Hellweg, Hans Horn, 
Christian Huber, Uwe Huniar, Marco Kattannek, Andreas Kohn, Christoph 
Kolmel, Markus Kollwitz, Klaus May, Paola Nava, Christian Ochsenfeld, 
Holger Ohm, Holger Patzelt, Dmitrij Rappoport, Oliver Rubner, Ansgar 
Schafer, Uwe Schneider, Marek Sierka, Oliver Treutler, Barbara Unterreiner, 
Make von Arnim, Florian Weigend, Patrick Weis, Horst Weiss. Available 
online at: http://www.turbomole.com. 



Im-[Co m (corrin)]-Rib + for consistency with previous studies 
(Jaworska et al., 2007; Lodowski et al, 2009, 201 1). In the RibCbl 
model (see Figure 1), the Ado group was simplified to ribosyl 
(Rib) by omitting the purine ring, leaving only ribose in the upper 
corrin plane. The corrin ring was truncated by replacing the side 
chains with hydrogen atoms, and DBI with imidazole (Im). The 
resulting overall charge of the system was + 1 due to removal of 
the PO4 nucleotide loop from below the plane. Previous studies 
have demonstrated that such simplified structural models accu- 
rately reproduce the electronic and spectroscopic properties of the 
full B12 structure (Jaworska and Lodowski, 2003; Jensen and Ryde, 
2003, 2009; Stich et al, 2003, 2004; Jaworska et al, 2005, 2007; 
Kuta et al., 2006, 2009; Liptak and Brunold, 2006; Kozlowski et al, 
2007, 2012; Rovira and Kozlowski, 2007; Galezowski et al, 2008; 
Lodowski et al, 2009, 2011). More inclusive structural models 
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FIGURE 2 | Absorption spectrum of RibCbl calculated in gas phase 
with BP86/TZVR 
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FIGURE 3 I Absorption spectrum of RibCbl calculated in water solution 
with BP86/TZVR 
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of coenzyme B12 will be employed in future studies utilizing 
multiple levels of theory. 

Optimized RibCbl structures correspond to equilibrium 
geometries where Co-Cr^ = 2.02 A and Co-Ni m = 2.21 A in 
the gas phase, and Co-CRib = 2.02 A and Co-Ni m = 2.19 A in 
water (PCM). Optimized axial bond lengths are in reasonable 
agreement with X-ray data of Randaccio et al. (Ouyang et al., 
2004) with corresponding experimental bond lengths of 2.04 A 
for Co-CAdo and 2.23 A for Co-Ndbi> respectively. Molecular 
orbital (MO) energies and fragment contributions of RibCbl are 
collected in Table S5 for the gas phase, and Table S6 for water 
(PCM). The relevant MO diagram can be found in supporting 
information as Figure S2. 

Frontier orbitals of RibCbl were determined by fragment 
analysis to be mainly corrin jt, jt*, and metal d orbitals. 
In particular, HOMO, HOMO-2, -4, -6, and LUMO+1, 
+2, +3 show major contributions from dc 0 and tt* while 
HOMO-3, -5, -8, -9, -10, and LUMO+3 have increased 
weights from the Rib group. On the other hand, LUMO and 
LUMO+4 are composed of nearly pure n orbitals of the cor- 
rin ring. Finally, HOMO-7 and LUMO+5 are both composed 
of excitations from the imidazole group. The change of envi- 
ronment from gas phase to water solution (Table S6) did 
not influence the character of the MOs significantly, differ- 
ences being attributed to reversion of energies between the 



HOMO-5 and -7 in the gas, and HOMO-3 and -4 in the PCM 
model. 

ELECTRONICALLY EXCITED STATES OF RibCbl 

Simulated absorption (Abs) spectra generated from TD- 
DFT/BP86/TVZP vertical excitation energies with transition 
dipole moments, using corresponding optimized RibCbl geome- 
tries in the gas phase and in water, are shown in Figures 2, 3 
respectively. Since hybrid B3LYP functional underestimates BDE 
as well as produce low-lying states as pure Tt colr i n —>■ it* ■ exci- 
tations, we employed BP86 throughout this study which has been 
proved to produce correct BDE for both MeCbl and AdoCbl 
(Kuta et al, 2006; Kozlowski et al, 2007, 2012; Galezowski et al, 
2008; Kornobis et al, 2011, 2013; Kumar et al, 2013). BP86 
functional was also frequently applied in QM/MM calculations 
for the enzyme-active cofactors (Kumar et al., 2012; Kumar 
and Kozlowski, 2013). In both the gas phase and solvent-based 
calculations, forty excited states were calculated to cover a range 
appropriate for electronic states potentially involved in coen- 
zyme B12 photolysis. The lowest ten singlet states are listed in 
Table 1, and the lowest ten triplet states are listed in Table 2 for 
the gas phase. In addition, the corresponding solvent-based cal- 
culations are listed in Table 3 for singlet states and Table 4 for 
triplet states. Complete lists of the forty excited states for each sys- 
tem can be found in the supporting information (Tables S1-S4). 



Table 1 | The lowest ten singlet states for RibCbl received from TDDFT/TZVP gas phase calculations. 
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Table 2 | The lowest ten triplet states for RibCbl received from TDDFT/TZVP gas phase calculations. 
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146 


H-1 


->■ L+1 


it + d xz /d z z -> d XK -n + 11* 








6 


144^ 


148 


H - 


• L+3 


d KZ + 11 -> a*(d z2 ) + n 








6 


143 


147 


H-1 


-> L+2 


it + d xz /d z z -» d XK -n + 11* 








4 


144^ 


147 


H -: 


• L+2 


d KZ + 11 -> d XK -n + it* 








4 


142 


147 


H-2 


-> L+2 


d X z + n -> d XK -n + it* 



Table 3 | The lowest ten singlet states for RibCbl received from TDDFT/TZVP PCM (water) calculations. 



E(eV) X (nm) f Coeff. Character 



Si 


2.28 


543.4 


0.0150 


27 


142 


145 


H-2 


-> L 


d xz + it -s- 11* 










56 


143 


145 


H-1 ■ 


-> L 


it + d yz /d z 2 it* 










15 


144^ 


145 


H -> 


L 


d yz + it -s- 11* 


S2 


2.39 


518.9 


0.0353 


74 


144^ 


145 


H -> 


L 


d yz + it -s- 11* 


S3 


2.45 


506.4 


0.0597 


61 


142 


145 


H-2 - 


-> L 


d xz + it -s- 11* 










29 


143 -»• 


145 


H-1 - 


-> L 


it + d yz /d z 2 -*■ it* 


s 4 


2.62 


473.5 


0.0024 


97 


141 -»• 


145 


H-3 ■ 


-> L 


d x 2_ K 2 -s- 11* 


s 5 


2.86 


434.0 


0.0043 


62 


144^ 


146 


H -> 


L+1 


d yz + it ->■ dxy -n + it* 










16 


140 -»• 


145 


H-4 ■ 


-> L 


nRib/cRib -»■ it* 


Se 


2.90 


427.7 


0.0235 


62 


140 -+ 


145 


H-4 


-> L 


nRib/cRib ->■ it* 










23 


143 -»• 


146 


H-1 - 


-> L+1 


it + d yz /d z 2 -> d xy -n + 


s 7 


2.94 


421.8 


0.0186 


53 


143 -»• 


146 


H-1 - 


-> L+1 


it + d yz /d z 2 -> d xy -n + 










14 


140 -»• 


145 


H-4 ■ 


-> L 


nRib/cRib ->■ t* 










14 


144^ 


146 


H -> 


L+1 


d yz + it — > d X y — n + it* 


Ss 


3.02 


410.0 


0.0053 


59 


144^ 


148 


H -> 


L+3 


d yz + it -s- o*(d z 2) + n 










20 


144^ 


147 


H -> 


L+2 


d yz + it — > d X y — n + it* 


s 9 


3.09 


401.4 


0.0185 


55 


144^ 


147 


H -> 


L+2 


d yz + it -s- d xy -n + it* 










19 


144^ 


148 


H -> 


L+3 


dyz + It -*■ o*(d z 2) + n 


S10 


3.16 


392.7 


0.0112 


35 


139^ 


145 


H-5 ■ 


-> L 


31 Im ~* t* 










33 


143 -» 


148 


H-1 ■ 


-> L+3 


it + d yz /d z 2 -> cf*(d z2 ) + 










15 


143 -» 


147 


H-1 ■ 


-> L+2 


it + d yz /d z 2 -> d X y -n + 










9 


138^ 


145 


H-6- 


-> L 


it + d yz -»■ It* 
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Table 4 | The lowest ten triplet states for RibCbl received from TDDFT/TZVP PCM (water) calculations. 





E(eV) 


X. (nm) 


Coeff. 










Character 


T 
1 1 


1 ~7K 


/Uo.o 


A Q 


144 —S 


1 /I E 

• I4u 


H -s- L 


dy Z + JI — ► It* 


T 

*2 


1.94 


637.8 


47 


143 —j 


■ 145 


H-1 


-s- L 


7t + dy Z /d z 2 — >■ Jt 


T~ 
'3 


O 1Q 

z. i y 


DO/. I 


/l R 


1/10 

I4Z — ; 


1 /i r 

• 1 4u 


H-2 


-s- L 


d x z + It — >- It* 


T, 


0 OK 
Z.ZO 


K/1 Q A 


QQ 
■JO 


144 — ; 


1 /I K 
• 1 4D 


H -s- L+1 


Uy Z + ji — > a X y — n + Tt 








/ 


1/1/1 
144 — i 


1 /I "7 


H -s- L+2 


d yz + tt — s> d xy — n + JI* 


T 


Z.4U 


0 lo.o 


J4 


1/lO 
I4Z — i 


1 /l c 
• I40 


H-2 


— >• L-f- 1 


dxz + ji — >• d xy — n + ji* 








6 


142 —j 


- 147 


H-2 


. 1 _i_9 


d X z + Tt -> d xy — n + ji* 


T 


Z.uU 


/i op; c 


0"7 

Z / 


1/11 
14 I — ; 


1 /i c 

• I4u 


j — | 3 


) [_ 


0 X 2_y,2 — > JI 








1 5 


141 —. 


- 146 


u 0 
n-o 


-V 1 -1-1 

— > L-f- I 


d x 2_ K 2 -> d xy -n + JI* 








0 


1/11 

I4 I — ; 


1 /I "7 
• 14/ 


u 0 
n-o 


_v I -l_9 
— * L-f-Z 


A ■.A n 1 w * 

d x 2_ K 2 — » a X y — n + n 


T 


Z.uZ 


yl QO "7 


1 0 
I z 


1/1/1 

144 — i 


1 /I "7 
• 14/ 


u 

n — 


* 1 -1-9 


dy Z + it — >■ d X y — n + n* 








1 "7 
1 / 


1/1/1 
144 —. 


1 /l D 
• l4o 


H — i 


► L+3 


dyz + — ► cr*(d z 2) + n 








IT 
0 


1/11 

14 I — ; 


1 A C 
• I4u 


H-3 


L 


d x 2_j/2 — >■ 








/I 

4 


1/11 

1 4 I — 3 


1 A K 
• 1 4D 


H-3 


-> L+1 


d x 2_^2 — >■ a X y — n + it 


T 
'8 


0 R/l 

Z.D4 


/1QQ O 


1 R 
I O 


1/11 
1 4 I — 3 


1 /I R 
• 1 4D 


H-3 


-> L 


™x 2 —y 2 n 








I O 


1/11 
14 I — 5 


1 A R 


H-3 


-> L+1 


A ^ A n 1 w * 

d x 2_ K 2 — * d xy — n + It 








O 


1/11 
14 I — 5 


1 A ~7 
• 14/ 


H-3 


-* L+2 


A ^ A n 1 w * 

d x 2_ K 2 — * d xy — n + It 








A 

4 


1/1/1 
144 —3 


1 A ~7 
• 14/ 


H -H 


► L+2 


dyz + it — > d xy — n + it* 


T 9 


2.61 


475.6 


15 


143 -j 


• 146 


H-1 


-> L+1 


it + d yz /d z 2 -> d xy -n + it* 








10 


143 -) 


• 147 


H-1 


-* L+2 


it + d yz /d z 2 -> d xy -n + it* 








9 


143 -j 


• 148 


H-1 


-> L+3 


it + d yz /d z 2 -> a*(d z2 ) + n 








6 


142 -j 


• 147 


H-2 


-* L+2 


d xz + it -> d xy -n + it* 


T10 


2.71 


457.6 


21 


142 -) 


• 148 


H-2 


-* L+3 


d xz + it -> cr*(d z2 ) + n 








13 


142 -) 


• 147 


H-2 


-* L+2 


d xz + it -* d xy -n + it* 








7 


143 -) 


• 148 


H-1 


-* L+3 


it + d yz /d z 2 -> a*(d z2 ) + n 



Taking into account that light exposure used in experimental 
studies was no shorter than 400 nm to induce Co-C cleavage 
(Walker et al, 1998a,b; Shiang et al, 1999; Yoder et al, 2001; 
Cole et al, 2002; Sension et al., 2004, 2005a,b; Harris et al, 
2007), in further discussion we will focus on the low-energy 
region. 

Since the studied RibCbl model is a truncated form of 
full coenzyme B12, no direct comparison to the experimental 
spectrum can be made. However, the predicted spectra can 
be correlated to that of AdoCbl with good agreement (Shiang 
et al., 1999). Several calculated electronic transitions can be 
assigned to experimental spectral bands. Specifically, the S2 at 
2.36 eV and the S3 at 2.37 eV agree well with the experimen- 
tal transition at 2.35 eV (528 nm), the S$ at 2.87 eV agrees 
with the experimental transition at 2.90 eV (428 nm), and the 
S10 state at 3.10 eV can be assigned to absorption observed at 
3.07 eV (404 nm). The same trend applies to PCM geometries 
(Table 3). Namely, the S2, S3, S6, and S9 can all be corre- 
lated with experimental absorption bands for AdoCbl observed 
at 2.35 eV (528 nm), 2.90 eV (428 nm), and 3.07 eV (404 nm), 
respectively. Furthermore, the simulated RibCbl Abs spectra are 
also in reasonable agreement with the simulated AdoCbl spec- 
tra published by Andruniow et al. (2009) further supporting 
the ability of the truncated AdoCbl model, (i.e., RibCbl) to 
accurately reproduce the experimental spectroscopic features of 
AdoCbl. 



The absorption spectrum of RibCbl calculated in the gas 
phase at the BP86/TZVP level is shown in Figure 2, and assign- 
ments for the lowest forty singlet states are collected in Table SI. 
Based on these results it can be deduced that the a/p bands are 
located in the low energy region from 2.00 to 2.80 eV, covering 
first five calculated singlet transitions, the D/E band ranges from 
2.80 to 3.30 eV and is composed of multiple excitations, and the 
high-energy region, or so-called y band, has significant inten- 
sity covering several transitions from 3.30 to 4.40 eV. Based on 
the nature of the calculated electronic excited states (Table 1), the 
a/f5 bands involve primarily dco/it contributions, the D/E 

bands result from various electronic transitions with more con- 
tributions from the Rib group, and the same tendency is observed 
in the y band, where oscillator strength increases significantly. It 
should be noted that ribose increases its role in photoelectronic 
transitions in the high energy region, i.e., the D/E and y bands, 
respectively. 

Table 1 summarizes the ten lowest calculated excited singlet 
states of RibCbl in the gas phase. As shown in Table 1, the 
Si excitation present at 2.24 eV, carries only a small oscillator 
strength (0.008). The major contribution from the Si state is an 
excitation from HOMO-1 to LUMO, in which n + dxz/d z 2 jt* 
is dominant. The MO diagram (Figure S2), shows that Si excita- 
tion also involves a transition from a a orbital of axial bonding to 
the it* of the corrin ring, and is identified as SBLCT (Lodowski 
et al., 2009). The second noticeable contribution from the Si state 
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is the transition from HOMO-2 to LUMO, identified as Jt + 
dxz — > it*- Based on the MO diagram it is apparent that LUMO 
is a pure corrin tt* orbital, to which all the excitations in the a/p 
band generally take place. The next two states, the S2 and S3, are 
adjacent closely in energy. The most significant components of 
the S 2 and S 3 states are HOMO LUMO and HOMO-2 
LUMO, where HOMO is a mixture of it and d yz and HOMO-2 
refers to Tt and d^. Therefore, the d to ti* transition is domi- 
nant in this region, and should be described as MLCT transition 
(Cole et al, 2002; Jaworska et al, 2007; Lodowski et al, 2009, 
2011). 

The S4 state has even less oscillator strength intensity com- 
pared to that of Si, in which both HOMO-3 and HOMO-4 con- 
tribute almost evenly to excitation at 2.58 eV. HOMO-3 derives 
from ribose, according to MO diagram, while HOMO-4 is com- 
posed from the Co d x 2_ y 2 orbital. Coincidently, the S5 involves 
the same set of orbitals which contribute to the S4 state. The only 
difference, considering their specific characters, is that the S5 state 
mixes a minor contribution from the HOMO-1— »■ LUMO transi- 
tion. However, the larger oscillator strength of the S5 state makes 
it the strongest component of the a/p band. 

Electronic transitions in the D/E band are more complex than 
those in the ot/p 1 band, and the oscillator strengths are noticeably 
smaller. This region can be represented by electronic transitions 
from the S6 to S14 with three strong excitations occurring at 2.87 
(S6), 3.10 (S10), and 3.21 eV (S12) (Figure 2). The S6 and S10 have 
dco/t tt* character while the S12 shows a drjo/rc CT co-C 
nature. The other small intensity transitions in this region are of 
similar character. 

More than twenty singlet states can be associated with y band 
ranging from 3.31 (S15) to 4.23 eV (S40) (Table SI). Among 
them, the S25 and S27 have the largest transition dipole moment 
and are composed of nRib/ffRib + d^ld^ + Tt -Hs-d^— n + tt* 
transitions with a minor dc 0 /Tt Tt* contribution. The S29 
is dominated by dco/Jt it* transitions, while the S32 has 
dco/t — »■ it* m character. Near the end of the simulated spec- 
trum, another noticeable transition of nRib/ffRib + dxz/d z 2 + it — > 
cr*(d z 2) + n type is found at 4.14 eV corresponding to the S37 state 
(Table SI). 

Since triplet excitations have zero transition dipole moments, 
we will not take them into account for spectral analysis. The 
first three triplet states are energetically below the Si (2.24 eV), 
where the Ti occurs at 1.72 eV from an electronic transition 
of HOMO^LUMO, while T 2 and T 3 are transitions at 1.94 
and 2.15 eV, originating from HOMO-l^LUMO and HOMO- 
2— >-LUMO, respectively. All three electronic transitions can be 
described as dc 0 /t — > Tt* type (see Table 2 and Table S2). 

The Abs spectrum of RibCbl calculated in water solution 
using the same level of theory was plotted in Figure 3. The 
energy of the y band remains unchanged in comparison to gas 
phase calculation. However, the a/p" band is narrowed due to 
the blue shift of the S5 to the D/E region. The Si, S2, and S3 
are assigned to dco/n — Tt* transition. A significant contribu- 
tion to Si state is composed of a tt + dxz/d z 2 — »■ tt* excitation 
and, like in the gas phase, is characterized as SBLCT. The S4 
has pure d x 2_ y 2 — > tt* character. The D/E band covers states 
S5-Si 3 . The states S5 and S7 are identified as d/ir — ► TT*/d— n 



excitations with a small admixture of long range CT (LRCT) type 
transitions (nRib/oRib — > tt*) (Andruniow et al, 2009). The S6 
state at 2.90 eV is a mixture of a dominant LRCT, and minor 
d/ir — »■ 7T*/d — n character. In contrast, the S7 state is mainly com- 
posed of a dco/it - ► Tt* transition, with minor contributions 
from SBLCT. The other intense singlet transitions in this region 
are assigned to the dc 0 /Jt -*■ Tt* transition, which is consistent 
with singlet excited states observed in the gas phase calculation. 
The complete assignment of electronic transitions is presented in 
Table S3. In water solvated system, the T1-T4 states are below 
the Si state (Table 4), and all possess contributions from the 
dco/Tt - *■ Tt* transition. 

MECHANISM OF PHOTOLYSIS 

In order to simulate the photodissociation of AdoCbl, the Co- 
CRib bond in RibCbl was elongated to mimic the light induced 
bond scission. The length of the Co-CRib bond was systematically 
stretched with the increment of 0.05 A in a range from 1.90 to 
2.65 A. For each constrained distance, the ground state geometry 
was optimized at the BP86/TZVP level, followed by TD-DFT cal- 
culations for both singlet and triplet states. Figure 4 depicts the 
potential energy curves of a manifold of excited states computed 
as a function of Co-Cra, bond length. 




— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 
R(Co-C) [A] 

FIGURE 4 | Potential energy curves of the lowest-excited singlet (red) 
and triplet (blue) states of the RibCbl model complex along the Co-C 
bond stretch computed at TD-DFT/BP86/TZVR The triplet repulsive state 
is denoted as 3 cr -> a*. 



www.frontiersin.org 



February 2014 | Volume 1 | Article 41 | 7 



Liu et al. 



TD-DFT study of AdoCbl photodissociation 



The S1-S5 singlet states, which belong to ct/p 1 band, have non- 
repulsive character, and the gap between the Si and higher states 
gets larger when the bond length increases. States denoted as S2 
and S3 are parallel at short Co-CRn, distance but cross at ~2.50 
A. On the other hand the S4 and S5 keep similar energy but visi- 
bly separate from the first group of transitions (S1-S3). However, 
at ~2.30 A these two states show avoiding crossing, and the S4 
approaches the S2 at long length bond region. The excited states 
higher in energy (S6-S10), display multiple crossings at even short 
Co-C distances. According to the TD-DFT data, low-lying sin- 
glet states do not participate in photodissociation because there 



is no 'Oco- 



°Co-c) dominant. Most likely such a transition 



requires an energetically high excitation and none of singlet states 



with predominant 1 (oco- 



J Co-C 



) character are in a reason- 



able transition region. Additionally, it is expected that such states 
would not dissociate to the cob(II)alamin and ribosyl radicals, 
but rather to ionic fragments. However, the low-lying triplet 
states appear to have <xco-c — > °c 0 -c character in the achievable 
energy region. By connecting the energy points of 0co-c-featured 
triplet states at different dissociation distances, the 3 (o"c 0 -c — ► 
°Co-c) state w ' tn re P u l s i ve character can be located (black curve 
in Figure 4). At 2.35 A, this state becomes the lowest level in 
terms of energy, and displays an avoided crossing with Ti state, 



that typically accounts for dc 0 /n it* excitation. Together with 
an increasing Co-Cr^ distance, the energy of the triplet state 
drops dramatically above a Co-C bond distance of 2.40 A, indi- 
cating that dissociation may not be properly described with a 
single-determinant wavefunction approximation at long Co-CRib 
distances (Kumar et al, 2011a,b). It is also reasonable to pos- 
tulate that there is a conical intersection existing between the 
two states. Although the current study has limitations on the 
longer bond distances, it is expected that the photolysis hap- 
pens at much shorter bond lengths, where 3 (crc 0 -c ~~ > ff Co-c) 
excitation is described properly. Further multi-configuration 
wave functions would be required to obtain a more detailed 
picture. 

POTENTIAL ENERGY SURFACES OF LOW-LYING EXCITED STATES 

Since the photoexcitation of B12 derivatives involves primarily 
the structural changes of the axial ligand involved in N-Co-C 
bonding, this investigation was extended to the PESs computed 
as a function of both Co-CRib and Co-Ni m distances. The Si 
and S2 PESs were constructed from the equilibrium structure 
of the RibCbl (So state), and systematically stretching the two 
axial bonds with an increment of 0.05 A. Geometries at each grid 
point were optimized by BP86/TZVP, and the vertical excitations 
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FIGURE 5 | Potential energy surfaces for singlet ground state and two lowest singlet excited states of RibCbl together with their vertical projections 
plotted as a function of axial bond lengths (expressed in A) calculated in gas phase with BP86/TZVR 
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were calculated by TD-DFT at the BP86/TZVP level of theory 
(Figures 5, 6). 

In gas phase PES contours (Figure 5), the So shows an equi- 
librium geometry with the lowest energy at 2.02 A for Co- 
Ciub and 2.21 A for Co-Ni m , while for the Si PES, the mini- 
mum occurs at the longer Co-Ciub bond (~2.17 A) but shorter 
Co-Ni m (~2.00 A) distance. According to pump-probe time 
resolved experiments, an intermediate state of AdoCbl is pop- 
ulated through photoelectronic excitation, leading to a partial 
mixture of adenosyl radical and cob(II)alamin. It is therefore 
reasonable to postulate that such an intermediate state can be 
identified as the Si, which according to TD-DFT calculations 
has the Co-CRib bond elongated and Co-Ni m shortened due to 
an electron density shift from the metal to the equatorial bonds 
of the corrin ring. The PES plot of the S2 is nearly parallel to 
the So surface. Estimates based on PES projections (Figure 5) 
predict the energy minimum of the S2 surface is located at Co- 
CRib ~ 2.04 A, and Co-N Im % 2.23 A, ± 0.02 A, compared to 
the equilibrium geometry. This indicates no significant change 
of axial bond distances occur from the So to S2 transition. If 
this is true, the energy difference between the So and S2 states 
are mainly attributed the geometry changes of the corrin ring. 



At this point, the distortion of the corrin ring, mainly due to 
7T — > it* electronic transitions, increases the potential energy of 
the system. Since the geometry change in the Si state involves 
Co-CRib bond stretching, it is likely that the Si state is the most 
involved in the photolytic cleavage of Co-C, and thus the for- 
mation of the radical pair generated when the Co-C bond is 
stretched. 

The PES of RibCbl in water retains nearly all of the general 
features of the gas phase PES (Figure 6). However, although the 
equilibrium geometry of RibCbl in water has similar axial bond 
lengths to those obtained in the gas phase calculation (with the 
exception of a slightly reduced Co-N bond distance of 0.02 A), 
there remain subtle differences in energies between the So and 
Si, and So and S2 states. The energy gaps between So and Si, 
as well as So and S2, are larger for the solvated geometry than 
those of the gas phase. The vertical energy gaps between the min- 
imum point of solvated ground state and the lowest two excited 
states are ~2.30eV (S0-S1), and ~2.40eV (S0-S2) respectively, 
while the vertical energy gaps in the gas phase calculation are 
~2.20eV (S0-S1) and ~2.35eV (S0-S2). It can be reasoned that 
the larger vertical energy gaps from the ground state to the first 
two excited states in the PCM calculation, that the polar effect 




FIGURE 6 I Potential energy surfaces for singlet ground state and two lowest singlet excited states of RibCbl together with their vertical projections 
plotted as a function of axial bond lengths (expressed in A) calculated in water solution (PCM) with BP86/TZVR 
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of solvent tends to increase the energies required for RibCbl 
excitation. 

SUMMARY AND CONCLUSIONS 

The purpose of this study was to explore the mechanisms of Co-C 
bond scission in AdoCbl upon light exposure using the TD-DFT 
method. To accomplish this we used a simplified RibCbl model to 
mimic the molecular features of photoexcited AdoCbl. The opti- 
mized RibCbl geometry at the DFT/BP86/TVZP level of theory 
was very close to the X-ray structure of AdoCbl with respect to 
the Co-C and Co-N ax bond distances, within the deviation of 0.02 
A. Electronically excited states were calculated with TD-DFT, and 
both low-lying singlet and triplet states were analyzed. Calculated 
singlet states of RibCbl were used to generate electronic spec- 
tra that agree well with the experimental UV-Vis for AdoCbl, 
especially considering the characteristic low energy a/p band. In 
addition, PESs were generated as a function of both Co-CRib and 
Co-Ni m distances for singlet excited states using TD-DFT. The 
Si state was identified as the key state of RibCbl photoexcita- 
tion, with contributions from both MLCT and SBLCT. Similar 
to the analysis of low-lying excited states of MeCbl, a repulsive 
3 (°Co-c ~~ °Co-c) triplet state was found to facilitate photon 
induced dissociation of the Co-C bond in RibCbl. 
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